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ABSTRACT 

We study the effect of pulsar rotation on timing of binary pulsars, with particular emphasis on 
the double pulsar system J0737-3039. Special relativistic aberration due to the orbital motion of 
pulsar changes both the longitude and latitude of the emission direction with respect to the pulsar 
spin axis. The former gives rise to a shift of the arrival time of the pulse centroid (this is the 
conventional “longitudinal” aberration delay), the latter results in a distortion (contraction or dilation) 
of the pulse profile on the orbital timescale. In the framework of the rotating vector model of pulsar 
emission, the amplitude of pulse distortion depends inversely on the variation of polarization position 
angle across the pulse. For small angle between the pulsar magnetic and spin axes, as inferred for 
PSR J0737-3039A from polarimetric observations, the pulse distortion is significant (~ 1%) and the 
associated “latitudinal” aberration delay is much larger than the longitudinal one. We show that 
by monitoring the arrival time of separate pulse components as a function of pulsar orbital phase, 
the latitudinal aberration delay may be easily measured with the current timing precision. Such 
measurement would constrain the spin geometry of the system. The latitudinal delay can also be 
detected by monitoring system’s orbital parameters on the geodetic precession timescale. Because 
of the near edge-on orbital orientation of the PSR J0737-3039 system, general relativistic bending 
of pulsar A’s radio beam near its superior conjunction also introduces spin-dependent time delays 
of similar order of magnitude as the aberration delays. In addition, light bending splits the pulse 
profile into two variable components, corresponding to two gravitationally lensed images of the source. 
Detection of lensing effects is challenging, but may be possible with existing technology. 

Subject headings: pulsars: general — stars: neutron - pulsars: individual (J0737-3039A, J0737- 
3039B) — gravitational lensing — binaries: general 


1. INTRODUCTION 

The recently discovered double pulsar system J0737- 
3039 (Burgay et al. 2003; Lyne et al. 2004) presents an 
unprecedented natural laboratory for testing our under¬ 
standing of general relativity and pulsar magnetosphere 
physics. One of the unique features of this system is 
the almost edge-on orientation of its orbital plane with 
respect to our line of sight. Currently the two best 
constraints on the system’s inclination come from the 
Shapiro delay (i = 88.7° ± 0.9°; Lyne et al. 2004; Ran¬ 
som et al. 2004) and from scintillation measurements 
(|z — 90°| = 0.29°±0.14°, Coles et al. 2004). Among other 
effects this high inclination leads to periodic eclipses of 
the millisecond pulsar (pulsar A) by plasma in the mag¬ 
netosphere of the normal pulsar (pulsar B) (see Lyne 
et al. 2004; Kaspi et al. 2004; McLaughlin et al. 2004; 
Arons et al. 2004; Rafikov & Goldreich 2004; Lyutikov 
& Thompson 2005). 

The edge-on orientation of the system also makes pos¬ 
sible gravitational light bending of the pulsar signal when 
the pulsar passes through its superior conjunction with 
respect to the companion. In Lai & Rafikov (2005), we 
studied the effect of gravitational lensing on the pulse 
intensity and the impact of light bending on the corre¬ 
lated scintillation measurements of both pulsars. We also 
evaluated the effects of light bending on the pulse arrival 
time: geometric time delay and the modified Shapiro de- 

1 IAS, Einstein Dr., Princeton, NJ 08540; rrr@ias.edu 

2 Department of Astronomy, Cornell University, Ithaca, NY 
14853; dong@astro.cornell.edu 


lay - both of these timing contributions are independent 
of the pulsar spin. 

Since pulsar signals are due to the beamed emission of 
a rotating neutron star (as opposed to radial pulsation 
of the star), additional spin-dependent time delays arise 
due to special relativistic aberration (Smarr & Blandford 
1976; Damour & Deruelle 1986) and general relativis¬ 
tic light bending (cf. Schneider 1990). The goal of this 
paper is to evaluate these spin-dependent delays in the 
PSR J0737-3039 system, and assess their detectability 
and potential for constraining the parameters of the sys¬ 
tem. We also examine the distortion of pulse profile due 
to aberration and lensing, and discuss an alternative way 
of analyzing pulse arrival time and introduce the concept 
of “latitudinal delay” that focuses on specific features in 
the pulse profile. For PSR J0737-3039, a measurement 
of this latitudinal delay may lead to useful constraints on 
the geometry of the system. 

In (21 we calculate the spin-dependent time delays due 
to aberration and lensing. In §3, we examine the pulse 
shape distortion/variation and calculate the latitudinal 
delays. The prospects of detecting these effects are dis¬ 
cussed in §4. 

Although the formulae derived in this paper are gen¬ 
eral for any binary pulsar system, we shall apply them 
to PSR J0737-3039, concentrating on the effects of ro¬ 
tational delays on the timing of the millisecond pulsar 
A, with pulsar B playing the role of the lensing com¬ 
panion. We adopt the following parameters for the sys¬ 
tem: M p = Ma = 1.337-Mq, M c = Mb = 1.25 M@, 
the spin period of pulsar A Pa = 22.7 ms, the or- 
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bital period Pb = 0.1023 d, the orbital semimajor axis 
a = 8.784 x 10 10 cm, eccentricity e = 0.0878, longitude 
of periastron u> = 73.8° (as of 2004), and the orbital in¬ 
clination angle i in the range between 90.14° and 90.56°. 


2. SPIN-DEPENDENT TIME DELAYS 

Because of special relativistic light aberration, the di¬ 
rection of emission in the comoving frame of the pulsar, 
N (|N| = 1), differs from that in the observer’s frame, n 
(|n| = 1), by the vector 

(An) A = n x (n x (3 p ) ~ n 0 x (n 0 x /3 p ), (1) 

where v p = c/3 p is the pulsar velocity relative to the 
binary barycenter 3 , and no is the unit vector from the 
binary system to the observer. Because of general rela¬ 
tivistic light bending, the direction of emission n at the 
pulsar position differs from no- This lensing effect is im¬ 
portant only around the superior conjunction of the pul¬ 
sar, giving rise to potentially measurable time delay pro¬ 
vided that the minimum projected separation between 
the pulsar and its companion is comparable to the com¬ 
panion’s Einstein radius (Lai & Rafikov 2005). The de¬ 
flection vector due to light bending (An) l = n — no is 
given by 4 


(An) L = 



R 

0 || 


( AR ± \ n 0 x (r x n 0 ) 

V R ) a|| ' ( } 


Here r = r p — r c is the position vector of the pulsar rela¬ 
tive to its companion, R = no x (r x no) is the projection 
of r in the sky plane, R = |R| = r(l — sin 2 i sin 2 ip) 1 ! 2 (ip 
is the true anomaly measured from the ascending node 
of the pulsar), r = a(l — e 2 )/(l + ecoscp) is the dis¬ 
tance between pulsar and its companion (<p = ip — u> is 
the orbital true anomaly measured from periastron), and 
0 || = asin*(l — e 2 )/(l + esinw) is this distance at the 
conjunction projected along our line of sight (a, e and 
u> are the orbital semimajor axis, eccentricity, and lon¬ 
gitude of periastron, respectively). Gravitational lensing 
gives rise to two pulsar images (specified by the subscript 
“±”); A R± is the displacement (in the sky plane) of the 
image position, R± = (A R±/R + 1)R, relative to the 
fiducial position, R, and is given by 


AR± = l 



(3) 


with the Einstein radius 

Re = (2R g a||) 1/2 ~ (2 R g a) 1/2 = 2550 km. (4) 

Here R g = 2 GM c /c 2 = 3.69 km and the numerical value 
is for J0737-3039A with pulsar B playing role of the com¬ 
panion. Combining eqs. 0 and 0 , the direction of 
emission in the pulsar’s comoving frame is N = no + AN, 
with AN = (An)^ + (An)^. The deflection angle of the 


3 The constant velocity of the barycenter itself does not affect 
timing measurement. 

4 Analogous to the classical work of Damour & Taylor (1992), 

our calculations of the lensing effects are performed in the static 
limit, which neglects the change in the position of the companion 
during the time it takes the photon to cross the binary. Including 
this propagation effect in the spirit of Kopeikin & Schafer (1999) 
would not qualitatively change our conclusions and we omit it from 
the present consideration for the purposes of clarity. 


emission vector due to light bending, Ai?/a, has a max- 
imum value (for the “plus” image) Re/ a ~ v/c (this 
occurs at the orbital conjunction and i = 90°), which is 
of the same order as the aberration effect. This similarity 
of the magnitudes of two effects holds only in the binary 
systems which are so highly inclined with respect to our 
line of sight that R ~ Re at conjunction (in more face-on 
systems with cos i ~ 1 the gravitational light deflection 
is a {v/c) 2 effect). 

Variation of the emission direction N causes a shift in 
the equatorial longitude 4> of N in the corotating frame 
of pulsar (counted in the direction of rotation) given by 


AN_(sp^no) 

|s p x no | 2 ’ 


(5) 


where s p is the unit vector along the pulsar spin axis. 
For an arbitrary emission pattern rigidly rotating around 
s p , any change in the equatorial longitude corresponds 
to a change of emission phase, leading to a time delay 
At = Adi/flp = (At )a + (A t)L (here is the angular 
frequency of the pulsar and the time delay is positive for 
signal arriving later), where 


Pp ■ ( S P X n o) 

(6) 

fip |s p x n 0 | 2 

AR± A r • (sp x n 0 ) 

(7) 

R ) Hpaylsp x n 0 | 2 

to aberration (Smarr 

& Bland- 


(At) A = 
(A t) L = 


ford 1975) and gravitational light bending (cf. Schneider 
1990) 5 , respectively. Note that this “longitudinal” time 
delay shifts the whole pulse profile uniformly without in¬ 
troducing any distortions to its shape (see §3). 

The pulsar spin axis s p is specified by two angles (see 
Fig. 1 of Damour & Taylor 1992; thereafter DT92): { is 
the angle between s p and the light-of-sight vector no, and 
rj is the angle between the ascending node of the orbit 
and the projection of s p on the sky plane. The aberration 
time delay @ is then given by 


(At )a = A (sin 'ifj + e sin uj) + B (cos ^ + e cos lo) , (8) 





n p cVl — e 2 sin £ 


( sin 7] 

1 cos i cos T] 


(9) 


where fib = ( GM t /a 3 ) 4 / 2 (M t = M p + M c is the total 
mass of the system) and a v = (M c /M t )a. These expres¬ 
sions coincide with the results of Damour & Deruelle 
(1986) and DT92 (expressed in terms of a complimen¬ 
tary angle A = tt — { instead of {). Similarly, the delay 
associated with light bending is 




r A sin rj cos ip — cos i cos r/ sin ip 
a\\J Sip sin £ 

/ A R± A (sin r) cos ip — cos i cos 77 sin ip) 
(ayflpy sin£(1 — sin 2 i sin 2 ip ) 1 / 2 


Note that the lensing effect is important only around 
the orbital conjunction (ip = 7t/2) for nearly edge-on 


5 This agrees with with expression of Schneider (1990) if we 
replace aii by r ■ no = r sin 7 sin p — this would lead to divergence 
at ip = 0,7r, which is naturally avoided in our case. 
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systems. Let i = ir/2 + Ai and ip = n/2 + Aip, with 
|A*|, |At/j| < 1 , then we have 

A R± \ A i cos r] — Aip sin 77 
a \\^p ) sin / [(Ai ) 2 + (Aip) 2 ] 1 ^ 2 
with AR± given by eq. © and R ~ 
0 || [(Ai ) 2 + (Aip ) 2 ] 1 ^ 2 . Both eqs. 0 and ITUl) dif¬ 
fer from the result of Doroshenko & Kopeikin (1995) 
(also quoted in Wex & Kopeikin 1999) since these 
authors used the equation for the light ray trajectory 
from Klioner & Kopeikin (1992) obtained under the 
assumption of R Re ■ As a result, in the most inter¬ 
esting case when R ~ Re (and bending delay attains 
its maximum value) the expression for (Af)^ obtained 
by Doroshenko & Kopeikin (1995) is not applicable and 
cannot be used for interpreting timing measurements 
6 . Also note that some previous treatments of bending 
delay (e.g. Goicoecliea et al. 1995) were restricted to 
circular and/or purely edge-on orbits. 

Neglecting the dependence on various angles, the aber¬ 
ration time delay is of order I2ba p /(cf2 p ) = 3.6/xs, and 
the rotational lensing delay at the conjunction is of or¬ 
der Ai?-t/(a||flp) = 10.5(A R±/Re) ^s. Not surpris¬ 
ingly, the two delays are of the same order of magni¬ 
tude at the conjunction, since flf,a p ~ v p /c ~ v/(2 c) and 

Re/ a ~ \[2vIc. 

Damour and Deruelle (1986) recognized that the aber¬ 
ration delay is degenerate with the Roerner delay, and 
therefore cannot be directly measured (but see §4). The 
time delays associated with lensing do not suffer from 
this degeneracy. Figures ^and | 2 ] show these delays (the 
combined geometric and Shapiro delays as calculated in 
Lai & Rafikov (2005), and the rotational lensing delay 
given by eq. m or EH) together with the amplification 
factor for both images, 

iR q. 0 1 

A± = W?+4 ± 2’ “=«/«£- ( 12 » 

In evaluating (A R)l, we adopt / = 50° as suggested by 
Demorest et al. (2004) from polarization measurement 
(although this angle is currently unconstrained) and use 
77 = 45° as an example in Figures Q]and [3 

One can see from Figure ^ that near conjunction the 
lensing delay associated with positive image reaches ss 
6 /xs for i = 90.28°. In binary systems with smaller incli¬ 
nation angles the timing signal due to lensing is weaker. 
One can show that if the minimum projected separation 
of pulsar from its companion a|| | Ai| Re then for e <C 1 
the maximum amplitude of (A i)^ is 

max|(At) L | » (1 + 1 cos 7 ? |), (13) 

21 Az | S sin( 

reached at Aip ss Ai(cosi] =p 1 )/sin 77 , where the up¬ 
per sign corresponds to cos 77 > 0 , while the lower one 
to cos 77 < 0. Note that increasing |Ai| diminishes 
max|(At)i| no faster than |Ax| _1 . The orbital phase 
at which this maximum occurs shifts from conjunction 
in proportion to |A*| reducing magnetospheric absorp¬ 
tion effects in systems with larger |A*|. In particular, if 



Longitude (degree) 

Fig. 1. The amplification (top panel), the combined geomet¬ 
ric and gravitational (Shapiro) delay (the second panel), the rota¬ 
tional lensing delay (the third panel), and the latitudinal lensing 
delay (the bottom panel) of the dominant (“plus”) image of the 
pulsar A signal as a function of the orbital phase. The longitude 
is measured from the superior conjunction of pulsar A (when A 
is exactly behind pulsar B). In each panel, the inclination angles 
are i = 90.56° (solid line), 90.28° (long-dashed line) and 90.14° 
(short-dashed line). For the third and bottom panels, the angles 
/ = 50°, 77 = 45° and tanxo = 0.08 are used. 


the real inclination angle of J0737-3039 is closer to that 
measured with the Shapiro delay (i = 88.7°) rather than 
the one inferred using the scintillation technique, then 
max|(At)i| amounts to ~ 1.7 fi s (which is at present un¬ 
measurable) reached « 13 s after conjunction, which is 
almost at the end of eclipse of pulsar A. 

The magnitude of (A t)^ scales as (sin </) -1 and its 
shape depends on 77 as (Ai)cos ?7 — (At/)) sin 77 . As a 
result, unlike the geometric delay and Shapiro delay, 
the rotational lensing delay is generally asymmetric with 
respect to orbital conjunction (the exception is when 
77 = 0, 180°). Note that although the delays associated 
with the “minus” image can become very large ( 100 ’s of 
/xs) as a result of its large deflection (|Ai?„/a|| |), this 
image is usually too demagnified to be of interest (see 
13T31 . 

3. PULSE PROFILE VARIATION AND “LATITUDINAL TIME 
DELAY” 

3.1. Pulse Profile Variation 

Variation of the pulse emission direction N also results 
in the change of colatitude of emission vector, £. Since 
cos/' = s p ■ N and sin/ = |s p x N|, we find 

AC = = ( A 0A + (A C)L, (14) 

|Sp X no I 

where 


6 By contrast, our expression is valid for all R as long as R <C a\\ ■ 
For R a||, the lensing effect is completely negligible. 


(A/)a 


-(n 0 x f3 p ) ■ 


(s P x n 0 ) 
|s p x n Q | 
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Fig. 2.— Same as Figure [Tl except for the subdominant (“mi¬ 
nus”) image. 


c\/l — e 2 


[cos i sin 77 (cos 4> + e cos lo) 


(A C)l = 


— cos 77 (sin t/> + e sin o>)], 
AR ± (n 0 x r) _ (s p x n 0 ) 
a\\R \s p x n 0 | 


(15) 



7* 

— (cos 77 cos t/j + cos i sin 77 sin t/>) (16) 
a \\ 


are the contributions from aberration and lensing, re¬ 
spectively. 

Unlike the longitudinal shift of the emission vector, 
the shift of colatitude affects the different components of 
pulse profile differently, thus changing the pulse profile 
periodically on the orbital timescale. In the following we 
consider for illustrative purposes a specific radio emission 
pattern based on the rotating vector model (RVM; Rad- 
hakrishnan & Cooke 1969) widely used for interpreting 
pulsar pulse profile and polarization data. In this model, 
the radio emission pattern consists of a set of cones which 
are circular in cross section as seen along the magnetic 
axis m (see Fig. 0 . As a result of pulsar rotation around 
s p , our line of sight n 0 passes through the emission cone 
with half opening angle p > |C — a\ (where a is the an¬ 
gle between s p and m) which causes two episodes of ra¬ 
dio emission at the pulse phases ±$o corresponding to 
the leading and trailing edges of the cone 7 . Emission 
patterns close to circular are supported by theoretical 
ideas of pulsar emission mechanisms and observations of 
slow pulsars (Rankin 1993), although their occurrence in 
millisecond pulsars is not so clear (Kramer et al. 1998; 
Weisberg & Taylor 2002). The observed pulse profile of 
PSR J0737-3039A is consistent with a circular emission 


' The radio intensities of the two emission episodes are gener¬ 
ally different because of the complicated azimuthal structure of the 
emitting region. 



Fig. 3.— Geometry of pulsar emission with respect to the spin 
axis s p , magnetic axis m, and our line of sight no- The latter cuts 
the celestial sphere and emission pattern (cone with opening angle 
p; the dotted circle denotes its intersection with celestial sphere) 
along the dot-dashed line in the course of pulsar rotation, giving 
rise to features of the pulse profile at ±$ 0 , shown schematically in 
the inset. 


pattern (Demorest et al. 2004; Manchester et al. 2005). 
The full width of the pulse is 24>o, with 

cos p — cos C cos a 


cos $ 0 = . . . 

sm C sm a 

The variation of £ then leads to a variation of 4> 0 : 

1 1 


A $ 0 = AC 


tan $0 tan C sin 4>o tan a 


(17) 


(18) 


It is easy to show that A4>o is given by (which is also 
clear from Fig. 0 


A $ 0 = - 


AC 

sin C tan xo ’ 


(19) 


where xo is the angle (on the celestial sphere) between 
the arc connecting no and s p and the arc connecting no 
and m at the edges of the pulse, i.e., it coincides with 
the position angle of linear polarization and is given by 
the usual expression (Komesaroff 1970) 


sin a sin 4 > 0 

tanxo =--T- (20) 

cos a sin £ — cos <Po sm a cos £ 

Equation ifTTJli shows that in the framework of RVM the 
change of pulse width depends inversely on the change 
of position angle of radio polarization across the pulse. 
It also shows that while AC is generally quite small 
(~ Vp/c ~ 10~ 3 due to aberration and ~ A R±/a ~ 
3 x 10 _ 3 A R±/R e due to lensing), A<h 0 can be signifi¬ 
cant when sinC tanxo is small, as in the case of the PSR 
J0737-3039 system (see M3 3.211 . 

Note that the variation of the pulse profile on the or¬ 
bital timescale due to the colatitudinal shift is not a 
unique property of RVM — this effect should exist for 
more complicated emission patterns as well. Interpreta¬ 
tion of a positive detection of this effect in terms of pulsar 
spin geometry would be less straightforward for intricate 
beam shapes, and the exact relationship between A$o 
and xo given by eq. would also break down, but 
they do not affect our basic conclusion about the exis¬ 
tence of such a phenomenon. 
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3.2. Latitudinal Time Delay 

Standard pulsar timing formulae (including those dis¬ 
cussed in § 2 ) measure the arrival time of the centroid 
of radio pulses. However, pulse profile often has multi¬ 
ple components; in the rotating vector model these com¬ 
ponents correspond to passages of n 0 through different 
emission cones (two for each cone). As discussed above, 
the variation of colatitude of emission direction due to 
aberration and lensing results in the change of the loca¬ 
tion (in phase) of a given pulse component relative to the 
pulse centroid. Thus if one times the arrival of a specific 
component of the pulse, additional delay associated with 
such shift of colatitude must be included; we call this 
“latitudinal delay”. 

In RVM, the latitudinal delay associated with the lead¬ 
ing edge of the pulse is 8 

(Af)( lat ) = A$o/H p = (A t)^ at) + (Af)£ at) , (21) 

where (A <)^ at1 and (Af)^ lat) are contributions due to 
aberration and lensing. Using eqs. m and (pi , we 
find the latitudinal aberration delay 

(At)^ at ' ) = C (sin ip + esiniu) + D (cos ip + ecosw)()22) 
where 


^ -1 -( _C0S ? ? 1 ,( 23 ) 

VlpcVT^ sm C tan Xo l cos * sm V J 

with tan X o( < f ) o) given by (12011 in the framework of RVM. 
Similarly, using eq. itRil . we find the latitudinal lensing 
delay 



(A t) { l at) = ( 


/ A R± \ r cos rj cos ip + cos i sin 77 sin ip 




\ R J ci|| sin £ tan X o 

( A R± \ A i sin 77 + Aip cos 77 


\a||Hp J sin ( tan X o [(Ai ) 2 + (At/;) 2 ] 1 / 2 ’ 


(24) 


where the second equality comes from expansion around 
ip = 7 T/2 + At p and i = 7 r /2 + Ah 

Naturally, if one times the arrival of a specific com¬ 
ponent of pulse profile, then the total delay is the sum 
of the longitudinal contributions [eqs. JS|, itTTl l associ¬ 
ated with the shift of the pulse centroid and the lat¬ 
itudinal contributions [eqs. E3, HI] associated with 
pulse dilation/shrinkage. Thus, the combination of lon¬ 
gitudinal and latitudinal delays results in ( 1 ) a uniform 
shift of the whole pulse profile in time by {At )a + (A t)^ 
and ( 2 ) a nonuniform profile contraction/dilation, with 
each pulse component being shifted in phase by A4>o = 
n p [(Af)^ at ' ) + (Af)^ at) ]. Pulse components closer to the 
pulse centroid (i.e. those with smaller <4> 0 ) are relatively 
more susceptible to contraction/dilation than the more 
distant components. 

To evaluate (Af)^ lat \ a knowledge of (/, \o and 77 is 
necessary. Analysis of polarization data of PSR J0737- 
3039A constrains the angle between the spin axis and 
the magnetic axis to be a k 4°, while leaving £ un¬ 
constrained (Demorest et al. 2004). The pulse profile of 
pulsar A exhibits remarkable similarity between its out¬ 
ermost leading and trailing edges (Demorest et al. 2004; 


8 For the trailing edge, the delay is —A$o/ttp- Hereafter our 
formulae refer to the leading edge of the pulse. 



b 1 

*C/3 

G 

<D 

G 



Ad> 0 



AO„ 


<J> 




m 



Fig. 4.— Variation of the pulse profile as a result of colatitudinal 
variation for two different cuts of our line of sight no through 
the circular emission pattern. The bottom panel shows the pro¬ 
jection onto the pulsar’s celestial sphere along the magnetic axis: 
the thick solid arc depicts part of the emission cone, the thin ver¬ 
tical arcs represent projections of meridians, and the dotted lines 
denote the magnetic field lines, a and b denote the locus of two 
different lines of sight. This Figure illustrates how the proximity 
of our line of sight to the edge of the emission cone affects the total 
pulse width. 


Manchester et al. 2005) at 4>o ss 115°, which suggests 
that all profile components come from the same mag¬ 
netic pole and that the two prominent spikes in the pulse 
profile can be associated with leading and trailing edges 
of the same emission cone. Based on the opening an¬ 
gle of other millisecond pulsars and £ ~ p, Demorest et 
al. (2004) suggests </ ~ 50°, which we adopt in our es¬ 
timate. Equation | 20 |) then yields tan Xo « 0.08. Such 
small value of tan X o results from the fact that in the 
Demorest et al. (2005) geometry our line of sight passes 
through the emission cone of pulsar A very close to its 
edge, see their Figure 3 for illustration. 

The bottom panels of Figures Q] and 0 depict the lens¬ 
ing latitudinal delay (Af)^* at ^ for both images. Clearly, 
because of the small value of tan X o for the PSR J0737 
system, (Af )^* at1 is much larger than (A t)^ and reaches 
« 79 /rs for i = 90.28°. In binary systems with smaller 
inclination angles satisfying a|||Ai| Re the maximum 
amplitude of (Af)^ reached at A ip « A*(l — sin 77)/ cos 77 
is given by 


max|(Af)^ at) | 


{Re/cl ||) 2 

21 A i | D p sin (/tan \o 


(1 + sin 77), 


(25) 


for e <?C 1. Scaling of max|(At)^ lat ' ) | with |Ai| is anal¬ 
ogous to that of max|(At)i|. For i = 88.7° and spin- 
magnetic orientation of pulsar A adopted here one finds 
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max|(At)^ lat) | « 22 /zs (which is likely still measurable) 
reached ft 13 s before conjunction, which is at the ingres 
of eclipse of pulsar A. 

Another constraint on the pulsar A spin orienta¬ 
tion (independent from the polarization measurements) 
comes from the observed lack of profile variation of A 
which should occur in general because of A’s geodetic 
precession (Manchester et al. 2005). Stability of A’s pulse 
profile suggests that A may currently be at a special pre- 
cessional phase with its spin axis almost coplanar with 
both the orbital angular momentum axis and our line of 
sight no, implying g ~ 90°. Assuming this Manchester et 
al. (2005) have determined the region of a-( phase space 
compatible with A’s pulse profile stability. A wide range 
of a is allowed, including the small-a solution preferred 
by Demorest et al. (2004). In Table |U we list the values 
of timing coefficients A , B , C, D , together with the max¬ 
imum amplitudes of longitudinal and latitudinal lensing 
delays (for i = 90.28°) for several spin orientations of A 
suggested by Demorest et al. (2004) and Manchester et 
al. (2005), assuming for illustrative purposes that in the 
latter case ij = 85°. Clearly, for small a (hence small 
tan\o), the latitudinal delays are much larger than the 
longitudinal ones, although the latitudinal aberration de¬ 
lay is reduced for ij ~ 90° (while the magnitude of the 
latitudinal lensing delay is largely unaffected). For illus¬ 
trative purposes, in the rest of the paper we will focus 
on the small-a solution (Demorest et al. 2004). 

3.3. Effect of Double Images 

In addition to the pulse profile change due to the 
change of emission colatitude (see &I 3. 111 . gravitational 
lensing affects pulse profile in another way. Near the 
orbital conjunction, light bending splits the source (pul¬ 
sar) into two radio images in the plane of the sky. The 
displacement of the image from the source is A R±, as 
given by eq. ®, and the image magnification A± given 
by eq. o. While aberration affects the arrival time 
of both images in the same way, lensing breaks this de¬ 
generacy and separates the arrival times of two weighted 
replicas of the unperturbed pulse profile. Assuming that 
in the absence of orbital motion and lensing the pulse 
profile is given by Io(4>), we find that the observed pulse 
profile at orbital phase ip, including the light bending 
effect, is 

/($; if) = A+ e~ T +I 0 ($ - D. p At + ) 

+A-e~ T ~I 0 ($ - n p At-), (26) 

where the factors e~ T± accounts for the absorption of 
pulsar signal by plasma in the magnetosphere of com¬ 
panion. The time delays A t± for both images include 
not only the rotational contributions [eqs. ©, 1101) . 12211 
and (EH)], but also the geometric delay and the (modi¬ 
fied) Shapiro delay (see Lai & Rafikov 2005). Note that 
A±, t- t, At± are all functions of the orbital phase ip. 

For illustrative purposes, we assume r + = r_, and 
show in Figure 0 the relative contributions of the two 
images in shaping the pulse profile /(<!>) near the orbital 
conjunction for i = 90.28°. In practice, one could adopt a 
simple procedure when analyzing arrival times near con¬ 
junction: since the “+” image dominates the lightcurve, 
one would naturally assume that the timing residual due 
to lensing is primarily given by At + and the negative 



Fig. 5.— (a) Pulse profile of J0737-3039A far from conjunc¬ 
tion, 7o(^ > )- (t>) Contributions to the total pulse profile f ( T) (solid 
curve) at conjunction (ip = 90°) due to the positive (dotted curve) 
and negative lensed image (dashed curve) taking into account dif¬ 
ferent time delays and magnifications of the images. In the calcu¬ 
lation i = 90.28°, ( = 50°, u = 45°, and a = 4° are assumed, (c) 
Difference between /(<1?) and /q(T), with the latter magnified by 
A+ + A and shifted by At+ in time at conjunction. The remain¬ 
ing residuals are due to the improperly subtracted negative image 
and are rather small. The flux units in the three panels are the 
same. 


image contributes to the pulse profile mainly through 
its magnification. Figure H: depicts the relative error 
incurred by this assumption: we plot the difference be¬ 
tween I 0 magnified by A+ + A- and shifted by At + in 
time and the real profile /($) which properly includes 
contributions from both images at conjunction. One can 
see that this simple procedure does perform quite well, 
and the flux residuals due to improper subtraction of the 
negative image are at most several percent. In reality, 
since the light of ” image passes closer to the com¬ 
panion star than the light of “+” image, and the mag- 
netospheric absorption depends strongly on the plasma 
density, we expect that near conjunction, t_ is much 
greater than r + . This would further decrease the con¬ 
tribution of the ” image to the distortion of the pulse 
profile. 

4. DETECTABILITY OF ROTATIONAL DELAYS AND PULSE 
DISTORTION 

The rotational time delays discussed in previous sec¬ 
tions could be observed in different ways. The aber¬ 
ration effect results in arrival time variation and pulse 
distortion throughout the whole cycle of the orbital mo¬ 
tion, while the lensing effect is important only in very 
edge-on systems near the superior conjunction of the 
pulsar. As discussed above, the pulse profile and po¬ 
larization of PSR J0737-3039A imply that our line of 
sight crosses the edge of the pulsar’s emission cone (De¬ 
morest et al. 2004), resulting in significant pulse distor¬ 
tion and the associated latitudinal delay. For exam- 
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pie, taking ( = 50°, tan xo = 0.08 and i = 90.28°, 
we find, for the aberration delays, A = —4.7 sin 77 ps, 
B = 0.023 cos 7 / ps, C = 59 cos 77 ps and D = 0.3 sin 77 77 s. 
The lensing delays are shown in Figures ^and[2] Aber¬ 
ration results in fractional pulse width variation of order 
A^o/^o = ^p(At)^ at V^o ~ flpC'/'ho ~ 10 _2 cos?7 in 
the course of the pulsar’s orbital motion. For a few sec¬ 
onds around the superior conjunction, light bending gives 
rise to additional profile distortion at a similar level. 

Motivated by our results, we propose to perform tim¬ 
ing of separate components of the pulse profile to detect 
pulse dilation/contraction and the associated latitudinal 
delays. Instead of fitting one template to the whole ob¬ 
served pulse profile to determine its overall timing dis¬ 
placement (or the arrival time of the pulse centroid), one 
can directly fit templates for each individual component 
of the profile. Assuming that some pairs of components 
represent the leading and trailing edges of the emission 
cone through our line of sight, one can measure the dif¬ 
ference between their times of arrival and correlate its 
variation with the orbital phase. This method presents 
a relatively easy way of measuring the latitudinal delays 
and extracting pulsar parameters from them. In the case 
of PSR J0737-3039A, periodic profile distortions must be 
very strong because of the small (~ 4°) magnetic incli¬ 
nation angle of pulsar A with respect to its spin (De- 
morest et al. 2004). The latitudinal delay is of order 
C = 59 cos 77 /lis, which should be detectable since RMS 
timing residuals of pulsar A are already at the level of 
27 p s (Lyne et al. 2004). Such detection would directly 
constrain xo, C and V- 

Detection of the standard “longitudinal” aberration 
delay (eqs. EHS]) is a more difficult task. Darnour & 
Deruelle (1986) were the first to recognize that this de¬ 
lay cannot be observed directly as a timing variation on 
the orbital time scale since it can be fully absorbed into 
the Roemer delay by rescaling/redefining the eccentric¬ 
ity, semimajor axis, epoch of periastron and two other 
post-Newtonian parameters. More precisely, because of 
the aberration delay, the observed values of the semima¬ 
jor axis a obs , eccentricity e obs , and the epoch of perias¬ 
tron Tq bs are related to the true values of these quantities 
gtrue^true^rptrue by the formulae 


d°b s 

e obs 




— 1 + eA, 


gtrue 

gtrue 


= 

cA 

« —3.3 x IIP 6 sin ? 7 , 

(27) 

a p sin i 

rjiobs 

J 0 

rptrue 
~ -*0 

- Vi ~ e 2 f2 V e Bi 


CB = 

cB 

1.6 x HP 8 cos 77 . 

(28) 


a p sin i 


where the numbers are for the J0737-3039 system. The 
degeneracy with Roemer delay implies that to detect 
aberration delay one must look for arrival time variation 
over longer timescales. As a result of geodetic preces¬ 
sion, both angles 77 and (, and thus a obs , e obs , Tg bs will 
change with amplitude ~ ca,£b, on the timescale of the 
precession period. In some systems the precession can be 
rather rapid (with period of tens of years) and, if ca, £b 
are not too small, this variation can be detected in long¬ 
term observations and the pulsar spin orientation can be 
inferred (Stairs et al. 2004). Such measurement would 


be complicated by the fact that variable aberration is 
not the only contribution to the evolution of a and e: 
they would also evolve because of the gravitational wave 
emission and possible time-dependent Doppler shift of 
the system, see DT92. However, these additional contri¬ 
butions should give rise to linear non-periodic trends in 
a and e making it possible to disentangle their effect in 
long-term timing observations from periodic (nonlinear) 
behavior of aberration time delay. 

In the case of the PSR J0737-3039 system, the preces¬ 
sion period of pulsar A is 75 yr, and the relative varia¬ 
tion of observed eccentricity and semimajor axis on this 
timescale is 3.3 x 10 ~ 6 sin 77 , while the observed epoch of 
periastron varies by 22 ps. The current timing accuracy 
of J0737-3039A (Stairs 2005, private communication) al¬ 
lows determination of e obs and a obs with relative accuracy 
of 1.7 x 10 -5 and 1.4 x 10 -6 , respectively, and the mea¬ 
surement of Tq bs has absolute accuracy of 35 ms. Thus, 
it would in principle be possible to detect variation of 
a °bs Q )Ut not of e obs or TQ bs ) caused by geodetic preces¬ 
sion and longitudinal aberration at a 2 cr level after tens 
of years of timing observation. 

Much better measurement can be made if one moni¬ 
tors the long-term drift of timing positions of individual 
pulse features instead of the whole profile. Then the am¬ 
plitude of drift 9 in a obs , e obs , Tg bs caused by geodetic 
precession is augmented by the much more significant 
latitudinal aberration. Since the latitudinal delay has 
the same dependence on the orbital elements as the lon¬ 
gitudinal one, its effect on variation of a obs , e obs , Ttf bs 
can be included by simply changing A —> A + C and 
B —> B + D in (1271) and (1281) . As a result, a obs and e obs 
determined from timing of individual pulse features vary 
on geodetic precession timescale with relative amplitude 
4.2 x 10~ 5 (cost7 — 0.08 sin 77 ). This should give a highly 
significant detection of geodetic variation of a obs (and on 
a shorter baseline than the longitudinal aberration vari¬ 
ation) and yield a measurement of e obs variation at a 2 cr 
level. The accompanying variation of bs at the level of 
0.4 ms would still escape detection. 

Time delays caused by gravitational light bending do 
not suffer from the degeneracy with Roemer delay (or 
any other timing contribution 10 ) because of the ir un ique 
time dependence (see eqs. m , m, and Figs. m - a s 
noted before, in general the rotational lensing delay is 
asymmetric with respect to the conjunction, which may 
facilitate the separation of this delay from the geometric 
delay and the (corrected) Shapiro delay (Lai & Rafikov 
2005). For the parameters adopted in Figure|T| the maxi¬ 
mum amplitude of (Af)^ for the positive image is ~ 5 ps 
(for i = 90.28° or the minimum R = 4000 km), while 

(Af)^ at ' reaches as high as 70 ps, both achieved within 
several seconds of conjunction. The latter delay would be 
detectable if one times the arrival of individual features 
of the pulse profile. 

Detecting the lensing-related time delays may be a 
challenging task. One factor which limits the accuracy 

9 In this approach the orbital parameters should be determined 
separately for each pulse component. 

10 The only degeneracy is with the frame-dragging delay caused 
by the companion’s spin (Wex Sz Kopeikin 1999; Rafikov Sz Lai 
2005) but in binaries consisting of two neutron stars the frame¬ 
dragging delay is always negligible. 







of pulsar timing arises from the intrinsic variation of the 
pulse profile from one pulse to the next. If one times the 
arrival of a pulse feature with width w, then the timing 
accuracy after combining N p pulses is a t ~ w/^/N p . For 
PSR J0737-3039A, the pulse has two peaks, each with 
width w ~ 4 ms. The lensing effect lasts about ~ 4 s, 
so that in each orbit there are N p ~ 200 pulses during 
the lensing event, and about half of them 11 may be used 
for timing purposes. This gives <r t ~ 0.4 ms for a sin¬ 
gle passage through conjunction, somewhat larger than 
the amplitude of (Af)^ lat ^ (for i = 90.29°). By adding 
up pulses from many orbits one will obtain better timing 

accuracy which may make the detection of (Af)^ lat ' pos¬ 
sible. This would set a robust constraint on the system’s 
inclination angle. 

Finally, all rotational effects (both aberration and lens¬ 
ing) that we described in application to J0737-3039A 
should also be relevant for J0737-3039B. Since the masses 
of the two pulsars are quite similar, the variation of the 
direction of emission vector AN is also similar for both 
pulsars, but the spin period of pulsar B, Pb = 2.8 s is 
significantly longer than Pa = 23 ms. As a result, the 
amplitudes of all spin-dependent delays for pulsar B are 
~ 10 2 higher than for pulsar A. However, the current 
timing accuracy of pulsar B, 2.7 ms (Lyne et al. 2004), is 
much poorer than that of pulsar A and is not expected to 
improve significantly in the future. The pulse profile of 
pulsar B is unstable and subject to variations through¬ 
out the orbit, presumably due to the influence of the 
intense wind/radiation from pulsar A (Ramachandran 
et al. 2004). Moreover, the spin and magnetic geome¬ 
try of pulsar B (not well constrained observationally at 
present) may not be very favorable for exhibition of large 
latitudinal effects 12 . 

5. DISCUSSION 

Pulse profile distortion associated with the latitudinal 
aberration delay is the most easily measurable among 
all rotational timing effects in J0737-3039A if the spin- 
magnetic orientation of this pulsar inferred by Demorest 
et al. (2004) is correct. This effect should be detectable 
with the current timing accuracy if the arrival times of 
individual pulse components (rather than the pulse cen¬ 
troid) are monitored. However, if this effect is not found 
in the data, we can suggest two possible reasons. First, 
it is conceivable that at the present time the spin ori¬ 
entation is such that cos?? <C 1, resulting in a rather 
small profile distortion by aberration (see eq. [ 231 )- This 
would additionally account for the remarkable stability 

11 Lensing effects manifest themselves only very near the con¬ 
junction (when R ~ Re) where absorption by magnetospheric 
plasma of pulsar B is important (Kaspi et al. 2004; Rafikov & 
Goldreich 2004). At the same time, McLaughlin et al. (2004) and 
Lyutikov &; Thompson (2005) have demonstrated that even very 
close to A’s superior conjunction radio signal of A can still pass 
through B’s magnetosphere almost unattenuated at some intervals 
of B’s rotation phase within which at least several rotational peri¬ 
ods of A can be detected, ft is this penetrating part of A’s radio 
flux that can be used to search for gravitational lensing effects in 
J0737-3039. 

12 A detailed modeling of the modulation of the pulsar A eclipse 
by pulsar B’s rotation yields the following parameters for pulsar 
B: the angle between the spin axis and magnetic axis ~ 75°, the 
inclination of the spin axis to the orbital angular momentum axis 
~ 60° (Lyutikov & Thompson 2005). 


of J0737-3039A pulse profile (Manchester et al. 2005). 
However, such a situation would not last since geode¬ 
tic precession of the pulsar spin axis changes rj on the 
timescale of years, so that several years from now latitu¬ 
dinal aberration should become observable. There is also 
a possibility that the spin of pulsar A is almost aligned 
with the orbital angular momentum 13 which, because of 
the edge-on orientation of the system, would naturally 
make cos rj <C 1. This requires a very large opening an¬ 
gle of the emission cone, p ~ 90°, which may be unlikely. 
Second, it is possible that the spin and magnetic orienta¬ 
tions inferred by Demorest et al. (2004) are incorrect, e.g. 
because of the failure of RVM. This would not be very 
surprising given that the polarization data are corrupted 
by quasi-orthogonal modes. Even the interpretation of 
the emission as coming from a single magnetic pole can 
be called into question 14 . If this were the case, the an¬ 
gles a and £ could be different from what we used in 
our estimates, possibly reducing the magnitude of pulse 
distortion. Thus, measuring the pulse variation in J0737- 
3039A on the orbital timescale can be used to probe the 
emission and spin geometry of pulsar A. 

Pulse shape variations caused by latitudinal shifts are 
very promising from the observational point of view: un¬ 
like the lensing effects, profile dilation due to orbital 
aberration does not require special orientation of orbital 
plane which makes it potentially observable in many bi¬ 
nary systems. As far as we are aware, no one has at¬ 
tempted to detect pulse profile expansion/contraction 
of the predicted form correlated with the phase of pul¬ 
sar orbital motion. The only possible exception is the 
study of PSR B1534+12 by Stairs et al. (2004) in which 
pulse shape variations were detected through fitting a lin¬ 
ear combination of two standard profiles to all observed 
profiles. Variation of coefficients of this linear combi¬ 
nation on short (orbital) and long (geodetic precession) 
timescales allowed Stairs et al. to set constraints on the 
pulsar spin geometry. This procedure seems to be less 
straightforward than a direct measurement of the varia¬ 
tion of the phase separation between different pulse com¬ 
ponents that we propose to perform, although it has an 
advantage of yielding pulsar spin orientation free from 
assumptions about the beam geometry. 

It may be possible to infer the PSR J0737-3039A spin 
orientation from rotational effects alone, without resort¬ 
ing to polarimetric observations. The pulse profile varia¬ 
tions on the orbital timescale and the long-term geodetic 
variations of orbital parameters may be enough to mea¬ 
sure A’s spin orientation with respect to its magnetic 
axis (angle a) and with respect to orbital angular mo¬ 
mentum of the binary. Detection of lensing effects would 
refine the accuracy of such measurement and addition¬ 
ally constrain the inclination of the system, and may also 
help verify the assumption of circular beam shape. 

Rotational delays may be measurable in other compact 
binary pulsar systems with significant v/c. In such sys¬ 
tems pulse profile variations on orbital timescales may be 
a telltale signature of latitudinal aberration. The effect 

13 The natal kick on pulsar B (e.g. Willems et al. 2004) would 
in general make s,; misalign with the orbital angular momentum 
axis except if the kick was directed along or close to the orbital 
plane. 

14 This would also affect constraints obtained by Manchester et 
al. (2005) which were derived assuming single pole configuration. 
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Table 1. Time delays for various spin geometries of J0737-3039A. 


a 

C 

V 

tan xo 

A, /is 

B, /is 

C , fj,s 

D , /is 

maxKAt)^!, /is 

max|(A£)^ at ) |, /is 

4° 

50° 

45° 

0.081 

-3.4 

0.016 

42.0 

0.20 

5.8 

72 

18.7° 

78.5° 

85° 

0.31 

-3.7 

0.0016 

1.1 

0.059 

3.1 

17 

5.7° 

83.7° 

85° 

0.091 

-3.7 

0.0016 

3.5 

0.20 

3.0 

56 

34.2° 

42.3° 

85° 

0.88 

-5.4 

0.0023 

0.54 

0.030 

4.5 

8.4 


is strongest for pulsars in which the polarization angle of 
radio signal does not exhibit significant variation across 
the pulse profile, implying small value of tanyo and large 
Ad’o (see eq. CHI). This could serve as a useful criterion 
for selecting potential targets to detect orbital pulse pro¬ 
file variations (geodetic variations will also be enhanced 
for such pulsars by latitudinal effects, see Timing of 
individual pulse components needed for such detection 
should be facilitated in pulsars with rich pulse profiles 
exhibiting many symmetric sharp features. Using this 
technique to measure pulsar orientation with high sig¬ 
nificance may require stacking together observations ob¬ 
tained in many orbits 15 , but, presumably, this would still 
require much shorter baseline than is needed for detec¬ 
tion of timing variations caused by geodetic precession. 

Our results on the lensing effects have been obtained 
in the static approximation, which neglects the change 
in the orbital configuration of the binary during the time 
it takes a radio photon emitted by pulsar A to cross the 
system. This is the approach adopted by Damour & 

15 Time interval during which such observations are taken must 
be small compared to the geodetic precession period since otherwise 
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